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NOON states (states of the form \N >a |0 >b +|0 >a \N >(, where a and b are single parti- 
cle states) have been used for predicting violations of local realism (Greenberger-Horne-Zeilinger 
violations) and are valuable in metrology for precision measurements of phase at the Heisenberg 
limit. We show theoretically how the use of two Fock state Bose-Einstein condensates as sources in 
a modified Mach-Zehnder interferometer can lead to the creation of the NOON state in which a and 
b refer to arms of the interferometer and iV is a subset of the total number of particles in the two 
condensates. The modification of the interferometer involves making "side" measurements of a few 
particles near the sources. These measurements put the remaining particles in a superposition of 
two phase states, which are converted into NOON states by a beam splitter if the phase states are 
orthogonal. When they are not orthogonal, a "feedforward" correction circuit is shown to convert 
them into proper form so a NOON results. We apply the NOON to the measurement of phase. Here 
the NOON experiment is equivalent to one in which a large molecule passes through two slits. The 
NOON components can be recombined in a final beam splitter to show interference. 



I. INTRODUCTION 

NOON states are interesting and useful they are "all-or-nothing" states, having the form 

\<i>) = ^[\N>,\0>b+\0>a\N>t] (1) 

where the subscripts a and b represent single particle states. Eq. ([T]) represents a superposition of all 
TV particles in state a and none in b, or none in b and all in a. Such states can have several important 
applications: 1) Since they are "Schrodinger cat" states of a system with particles, one might use them 
to demonstrate the quantum interference of macroscopically distinct objects f^. 2) They can be used to 
study violations of quantum realism in the GHZ contradictions l^i^. 3) They can be used to violate the 
standard quantum limit and attain the Heisenberg limit in metrology by providing extremely accurate 
measurements of phase 4) They have been proposed for use in quantum lithography [7]. NOON 
states have been made experimentally with N up to ten particles |8l-ll2|. A two-body NOON state can 
be constructed by allowing two bosons to impinge on either side of a 50-50 beam splitter. This is because 
the final state will be a sup erposition of two-particles on either side of the splitter according to the 
Hong-Ou-Mandel effect We have shown previously how one can use two Bose-Einstein condensate 
Fock states as sources for an interferometer that can produce NOON states, with a and b two arms of 
the interferometer [ij, [l^ . 

The essential idea of the interferometer is as follows. By drawing off a portion of our condensate 
particles into a pair of detectors Dl and D2, we measure phase, which puts the uncounted particles in a 
double phase state-a Schrodinger cat. If the phase difference of the two cat branches is tt, then passing 
the remaining particles through a beam splitter results in a NOON state in the two output arms of the 
interferometer. However, this result occurs only if the number of particles initially detected is ideal (equal 
numbers in Dl and D2). The phase difference can then be adjusted by use of a correction circuit, in 
which a "feedforward" method sets the transmission coefficient of a side detector D9 to an appropriate 
value The result is that a very good approximation to a pure NOON state can almost always be 
generated whatever the count of the initial detections at Dl and D2. 

We begin in Sec. |TT]by introducing the two-stage interferometer, with several parameters to be deter- 
mined. In Sec, mil we look first at the case that Dl and D2 register the same particle number, and show 
how a NOON state is obtained at the output by considing the distribution for the relative phase arising 
from the initial measurement. In Sec. IIVI we progress to the general case, for which Dl and D2 detect 
different particle numbers. To obtain a NOON state at the output, a condition is derived relating the 
transmission coefficient for D9, and the outcomes at Dl, D2 and D9. To address the probabilistic nature 
of the side detections, the average value for D9 is computed is Sec. |Vl providing a simple value for the 
transmission coefficient. In Sec. IVIl the scheme is evaluated by applying two NOON quality factors to 
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the states at the output. The efficiencies of the scheme without and with the correction circuit, tak- 
ing into account all possible measurement outcomes, are compared and explained in Sec. IVIIIl Finally, 
applications of the scheme to phase estimation and to demonstrating quantum interference are given in 



The interferometer to be used is shown in Fig. [T] Two Fock state sources of number Na and A^^ enter 
the interferometer. We will see that the side detectors 1 and 2, situated immediately after the sources, 
are a key element; by measuring mi and m2 particles in these detectors, the uncounted particles, in arms 
3 and 4, are put into phase states. If the phase relation is correct (equal numbers mi and m2 in detectors 
Dl and D2) then, when these remaining particles pass through the middle beam splitter, the result (for 
suitable value of ^) is a NOON state in arms 5 and 6. However, if toi ^ m2 then we will show that the 
transmission coefficient at detector D9 can be adjusted to a value that corrects the relative phase giving 
a NOON after the beam splitter at 7-8. 



FIG. 1: An interferometer, in which sources have numbers Na and A'^^, which are partially diverted to interfere at 
a beam splitter and detectors Dl and D2. The counts in these detectors are used to determine the transmission 
coefficient in the beam splitter leading to D9 correcting the phases of the beams in arms 5' and 6 so that an 
approximate NOON state emerges 7 and 8. Three phase shifters 8, ^, and must be properly adjusted. 

The annihilation operators at the detectors are found by tracing back from from a detector to each 
source. We have 
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II. INTERFEROMETER 
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where 




(3) 



and T and R = 1 — T are the transmission and reflection coefficients at the beam splitter leading to D9. 
We will immediately take 9 = 7r/2, which puts the phase states symmetrically around the zero angle as 
we will see. 
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III. CASE WITH NO CORRECTION CIRCUIT 

We first look at the uncorrected situation where we select only cases in which mi — ra^- We then 
want to show how a NOON state arises in arms 5 and 6. The amplitude for finding particle numbers 
{rni, 7712, ms, rrig} in those detectors is 



mi ,m2 ,m5 ,mfi 



^5 "6 "1 "2 



Put in the forms from Eq. ^ and expand the binomials to give 



(4) 



7712,7715,7716 



{p.} 



s(:;)(::)(::)(;":)'-"™-"<"'='""-<-""- 



X (^0 |aPi+P2+P5+P6g™i+™2+rn5+m6-pi-P2-P5-P6 I N^Np"^ 



(5) 



The second line is evaluated to -y/ NalNfj\Sp-^-^-p2+p^+pg^Na^m.i+m2+m5+me,N : where N = + Np. For 
accurate analysis we could replace one of the summation variables by use of the ^-function. But for 
physical analysis we replace the first by 



"Pl +P2 +P5 +P6 , -'Vq 



JiPl+P2+P5+P6-Na)<t> 



2tt 



(6) 



Putting this into Eq. ([5]) we find that each Oq. has been replaced by e**^ and each by 1. The result 
of redoing the sum is then 



C, 



7711 ,TTT-2 ,fTT'5 ,^6 



2^ y/mi\m2\m^lmQ\ J-tt ^tt 



(7) 



where 



Factor out e"^/^, to give 



where M — mi + m2 and 



(9l2(0) = I cos 



(8) 
(9) 

(10) 



For arbitrary mi, m2, the function Q12 has peaks at ±0o — ±2 arctan( •\/m2 / mi ) . In the case toi = m2 
the plot of Q12 has peaks at ±7r/2 as shown in Fig. [21 

To test for the presence of the NOON state we approximate the peaks in Q12 of Fig. [2] by (5-functions: 



The result is then 



C, 



7711 -7712 ,7715 ,'>T^6 



Qi2{4>) - 5(0 - 0o) + i~ir'Si^ + 0o) 



(11) 



(12) 



When we take ^ = and 0o = 7i'/2 first line is proportional to 0™'^(— 2)™** requiring ms = and 
mg = N—M for non-zero contribution. The second line is proportional to 2™'so™6 requiring mg = 0. Thus 
we get a NOON state with a superposition of |mim2m5mg) = |M/2, M/2, 0,N — M) and |mim2m5mg) = 
\M/2, M /2, iV - M, 0). To make 0o = 7r/2 we must have mi = m2. 
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FIG. 2: Plot of Qi2{4>) of Eq. (1101) for mi — m2 = 15. Because we have taken 9 = 7r/2 and mi = m2 we find 
the two peaks at ±7r/2. For odd m2 we have one positive and one negative peak. With rn2 even both peaks are 
positive. 

To avoid a numerical integral, we use the (5-function generated in Eq. ([5]) to keep the probability in the 
form of a sum: 



P, 



E 



-»(p+9)(€+7r/2)(^_2)P+'' 



p!(mi — p)\q\{m2 — q)\r\{m^ — r)\ 



{Na-p-q- ry.{p + q + r + me- iV„)! 



(13) 



By adjusting the phase to ^ = we get the approximate NOON state as seen in Fig. [31 The NOON state 
is not perfect because of the finite width of the peaks in Qi2. As rrii and m2 increase the peaks narrow; 
in the limit in which they can be replaced by (5-functions, the output state becomes an ideal NOON. 
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FIG. 3: Plot of Pmi,m2"i5,»"6 ^q. (O versus ms (me = 60 — ms) for Na = Nfi — 30, jni = m2 = 15, ^ = 
(solid line). We also show mi = 18, m2 = 12 (dotted line) to illustrate a case where the phase states are not quite 
orthogonal. Discrete points are conneced by lines as a guide to the eye. 



Thus we see in this case that we need mi to be only approximately equal to m2 to get a good 
approximate NOON state. 



IV. THE CORRECTION CIRCUIT 



Consider now the complete circuit of Fig. [TJ What do we have to do to produce a NOON state in D7 
and D8 in the case when mi ^ m2? The basic idea to answer this question was given in Ref. 14. The 
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amplitude for finding the detector counts {mi, m2, my, ms, mg} is 



,7712 ,m7,ms ,m9 — \ 



Jj'j ttg ttg Ui-y Ui'2 



\/milm2\mT\mslmg\ 



Pi J \P2 J \P7 J VPs 



i \ / \ f^-.m2-p2+ms-ps 



E 

^Spi+P2+P7+P8+Pa,Na^mi+m2+mr+m8+ma,N 



(14) 



Replacing the 5-function by an integral as above gives us 



^1 ^1^2 ,^^7,*Tl8 ,mQ 



2^V2'"'^'"Vmi!m2!m7!TO8!m9! J-^ 



■^-^129(0) 



where 77 is an unimportant phase and 

Ri29i(t>) = (e"^ + l)™i(e''^ - 1)™^ (ze'«e"^ + l)"' 



(15) 



(16) 



We see immediately that if we take ie^^ = 1 then the 9-terni has the same binary form as the 1-term and 
mi and uig will simply add. With ^ = — 7r/2 we have 



R 



129 



2mi+m2+ni9 jm2gi(™i+m2+m9)f | cOS — I 

V 2; 



mi+mg / j_ \ 7712 

sin ■ 



(17) 



which has peaks at ±0o — ±2 arctan(-\/'Ti2/(mi + mg)). 

However, for arbitrary (j> the second line of Eq. ([15]) becomes (for e'^ = —i), 



V 1 cos I sm — 

2 2 



(-iue*'^ + w)™'(ive^'^+M)™' = (-i)™^ VTe<(e*'^ + l)-(e^'^-l) VTe<(e*^ + 1) + (e''"^ - 1) 

— ^__j'jm7gj(7ri7+7ra8)0/22'™7 + 77lS 

/ — it <t> (t> 
X V Te^ cos - + i sin - 

We can make the factors real if we take ^ = 7r/2, which gives 



(18) 



a 



mi ,Tn2 ,m7 ,m8nig 



!TO2!m7!r7ig!m9! •^-•t 



^e'(^^-^°^*/^Ol29(</>)Q8(0) 



T cos — — sin — 
2 2 



where = toi + 7712 + mg, and 



(3i29((/') = ( cos- 



7711+7719 



sm- 



(19) 
(20) 



Tcos- 



(21) 



Q129 peaks sharply at two angles analogous to Q12 of Eq. (jlOp and, for large mg, Qs peaks sharply at a 
different angle; the product also peaks sharply at an intermediate angle. The product Q129QS is plotted 
in Fig. [3] for some specific parameter values to show this. 

If m7 is small, we can assume the last m7 factor in Eq. (|18p (call it A(0)™^) is slowly varying compared 
to the other factor, which can be approximately represented as a ^-function at its maximum angle 
Our result for the probability will then have a factor A{(j)m)™''' ■ If we make A((/>,„) as small as possible, 
then the probability will be small for all m7 except for m7 = 0; that will give us a NOON state. 
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FIG. 4: Plot of Q129QS versus <j!> showing that it peaks sharply. The location of the peak is given by the cubic 
equation given in Eq. (|22p . Here mi = 22, 7712 = 8, mg = 14. 

We can find the maximum angle of Qi29Qa by taking the logarithmic derivative of the quantity. We 
take 77178 equal to the total number of particles entering D7 and D8, that is mjs = mj+mg, = N ~mi—m2. 
The result of this is the following cubic equation in X = tan ^ : 

{mi + mg)X^ + VT{nii + mg + m7s)X^ - (7712 + mjs)X - m2\/T = (22) 
Also we want A{(j>m) = or 

Vt = X (23) 



Combining the last two equations gives us a value for T: 



2(7771 + mg) + 77178 

Since 7^78 = N — mi — 7712 — mg we have 

T = N~{m,-m2 + m,) 

N + (toi — 7772 + mg) 

If 77ii > 7772 this value of T will be surely less than one. However, if 7712 > ttii then we will have to go 
back to Eq. (|16p and take ^ = +7r/2 to give an appropriate T value, which has ttii and 7712 interchanged. 

Unfortunately, our value of T depends on the value of 7779, whose probability distribution in turn 
depends on T. We consider how to solve this problem in the next section. 



V. FINDING mg 



Our expression for T contains 7719, which we should approximate in some way to get the best NOON 
state. We want to set the transmission coefficient to a value that depends on the count that just occurred 
in Dl and D2; of course, then the number going into D9 is probabilistic and is not precisely known. 
However we could hope to do well enough by replacing 7779 in Eq. (|25[) by (ttiq) so 



T = 



N - {mi - 77i2 + {mc,)) 

N + (7771 — 7772 + ("79)) 



(26) 



We then need to find (for fixed mi and 70,2) ("T-g), which itself depends on T. 

We proceed just as before to now look at the probability of finding m^i ,me,mg particles just before 
the last beam splitter in Fig. [1] with given input values of ttii and 7712. The amplitude for this is 



nil m2^"»5' „ni6 _m9 



a 



flg a, 



9 



\/mi\m2lm5ilmQlmg\ 



(27) 
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We now note that, with the phases chosen, all of the operators are of the form ±a^) and for mi > 1x12 
we have 



m,i,m2 ,7710 ,7719 



2^v 



(aa+aff)'"^+"'-^'+"°(aa-a/3) 



7712 4- "/Tie 



(28) 



Using this relation we are able to show the following rigorous relation for the average of mg over a series 
of measurements at fixed mi and 7712: 



(mg) = (1 - T) (iV - mi - m2 - (me)) 

Particle conservation requires 

(ms) + (mg) — N — mi - m2 
Further we can prove that an extremely good approximation is given by 

(ms) sa — (mg) 
m2 



(29) 
(30) 

(31) 



We derive these equations in Appendix A. The physically revelant solutions to these equations, valid for 
mi > m2 are 



[mg) 



1712 

mi 

mi — m2 
mi + m2 



(N — mi — m2) 



(32) 
(33) 



We can actually compute a rigorous numerical average for (mg) (Appendix B) to show that this formula 
is quite accurate. 

VI. THE CORRECTED NOON STATE 



Now let us return to the expression for the NOON probability and compute the distribution with the 
most probable value of mg now known. One possible formula uses the 5-function in Eq. (jl4p to eliminate 
pi and to find the following result for the probability: 



mi ,m2 ,m7 ,ms ,mg 



i?™''(mi + m9)!2m7!m8!A^Q!iV/3! 
4^2™7+m8TO^!mg! 



E 



(-1) 



P2 



P2-P7 



(mi +mg +P2+P7 +Ps- Nay-P2Km-2 -P2)!p7!(™7 -P7)!P8!('^8 -PsV- 



(34) 



This formula is very general, but has many sums. Nevertheless we can use the angular expression given 
in Eq. (fT9|) . which is also exact and leads to fast and accurate numerical calculations. We have then the 
alternative formula 



P, 



R'^oNJNpl 



1^1 ,m7,m8 ,mg 



2™''+"^smilm2lmjlms\ 



vTcos ^2 2' 



27r 



/7f (P ■ V 

V J COS sm — 

2 2 



mi+mg 



(35) 



This probability is symmetric in exchange of mj and mg; to see this, after the interchange, simply set 
(j) — and the result is the same. We summarize the procedure: For a given mi, m2 we determine T 
and (mg) from Eqs. (|32p and ([55)1 . The value of mg used in the probability (other than in the T and R 
values themselves) should be at or near the most probable value (rounded to the nearest integer) since 
the probability distribution is fairly narrow (see Appendix B). 

In the first trial we pick the most probable mg value and then a less probable value of mg with use of 
the optimum T value. In the case shown in Fig. [S]e, Na = Nr = 35, mi = 22, m2 = 8, the value of the 
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FIG. 5: Pmi,m2mT,mgms versus 7717 wheii = — 7r/2, C = 'r/2 for: (Left) = = 35, mi = 22, 7112 = 8 , 
mg = 18 (the most probable value), T = 0.37, the value from Eq. (|26|) . Here the NOON qualities are qi — 0.97 
and 52 = 0.99. (Right) Same parameters except with mg = 14. We find qi = 0.91 and 92 = 0.98. 



transmission is T = 0.36 and the average value of (mg) — 18.2. When mg — 14, that is, a less probable 
value (the probabihty of getting this value relative to the most probable value is ~ 0.4), the distribution 
is only slightly less NOON-like. 

There is a simple NOON quality factor to test the approximate NOON state, namely 



91 — '2Pmi,m2,0,N~ 



-rni — 7712 ^77lg ,7779 



(36) 



that is, twice the value of the probability at my = 0; we would like qi to be as close to 1.0 as possible. The 
NOON quality for the mg = 18 case in Fig. [S]is 0.97; for mg — 14 it is 0.91. Even when the probability 
of having my = or myg = N — mi — m2 — mg is high, it is possible to imagine a peculiar situation where 
the rest of the particles might be situated near ■m7s/2, which would diminish the NOON quality of the 
state. A quality factor that takes this into account is 



92 = 4 



A^mr 



(37) 



'78 



where where A^my 



(my)^ is the variance of my over the distribution. In a perfect NOON state 



the variance is maximal and 172 ~ 1- The worst possible case might be when P{mj = — 1) = P{m7 = 
+ 1) = 0.5 with all others zero; this case has g2 = 0. A case in which all probabilities are equal 
has q2 = 0.33. The second quality number seems less sensitive to changes in the resulting NOON states, 
but compare them in the case shown in Fig. [HI While the first quality factor vanishes, the second one is 
greater than zero because there is some NOON-like separation in the two peaks. 



Prob 




FIG. 6: Probablity distribution for a poor quality attempted NOON. Here — 100, mi — 35, m2 = 5, T* = 1 
and mg = 0. The NOON quality factors are qi = and 92 = 0.34. 



Next we choose a wide range of values of mi , m2 to provide the thorough set of comparisons. The result 
given in Table HI which lists values of the optimal transmission coefficient at N = 140. The parameter 
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TABLE I: Quality factors for a wide range of input mi and m2 values, with the resulting transmission coefficients 
and mg averages for A'' — 140. 



mi 


m2 


m78 


(mg) 


T 


qi 


92 


45 


5 


18 


72 


0.11 


0.976 


0.990 


40 


10 


36 


54 


0.25 


0.968 


0.993 


35 


15 


54 


36 


0.43 


0.955 


0.993 


30 


20 


72 


19 


0.67 


0.932 


0.992 


25 


25 


90 





1.0 


0.883 


0.988 



77178 = + nis is the total number of particles involved in the NOON state. In each case the average 
value of 777g is used. 

It is likely, in any set of experimental runs, that a random assortment of values of 77ii, m2, and 7779 will 
be averaged over in making a NOON state. What percentage of the inputs will result in good NOON 
states? In Appendix C we compare the corrected and uncorrected efhciencies and show that the correction 
process is successful in producing good quality NOON states with high probability. 



VII. APPLICATIONS OF THE NOON STATE 



A. Metrology 



A key application of NOON states is to ultrasensitive sensors, with fundamental sources of noise 
reduced to the minimal level permitted by quantum mechanics [l^. Fig. [7Ii) shows a Mach-Zehnder 
interferometer set up for the measurement of a path-length difference X) a model that can applied to the 
detection of a variety of physical parameters. A wide variety of input states, measurement protocols, 
and decoherence models have been considered in the literature. In what follows, we will consider the 
usefulness of approximate NOON states generated by the interferometric method described previously. 
The phase estimation process can be considered as an additional interferometric stage at the output, as 
illustrated in Fig. [Tt^ii). 




detection 
observable 




FIG. 7: The application of our NOON-state generator to phase estimation is illustrated. The aim is determine 
the unknown path-length difference in a Mach-Zehnder interferometer with precision approaching the Heisenberg 
limit. In (i), we show the usual Mach-Zehnder setup for measurement of phase. In (ii), an additional stage is 
added to the NOON generator for phase measurements. The NOON is in arms 7 and 8. 



A typical phase estimation procedure is as follows. For each estimate of parameter with value 
X, the experiment is repeated t times. It is assumed that the (conditional) probability distributions 
P{k\x), for the possible detection outcomes k, are known from theory. Following a Bayesian approach, 
the posterior distribution P {x\kt ■ • • ki) is obtained from the prior distribution P (x) by the update 
rule P {x\kt ■ ■ ■ ki) cx P {kt\x) ■ ■ ■ P {ki\x) P ix) ■ ^ suitable estimator - such as the mean value for 
P{x\kf--ki) - is applied to obtain a value for x- For a so-called "global" estimation procedure no 
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prior information is assumed and P{x) oc 1/2tt; for a "local" estimation procedure, the aim is to track 
small changes of x [13 • It is for the latter case that NOON states are particularly useful. A common 
measure of the statistical information available for a local estimation procedure is the (classical) Fisher 

information defined as Id = J2k ^(^Ix) In [-P (^Ix)]) • Given a total of ly independent estimates 
(Xij ■ ■ ■ iX;^)> the Cramer-Rao bound places a lower bound on precision Ax (defined as the root mean- 
square error of the final estimate) as Ax > \l\fvtT^\. This bound may be assumed to be tight, provided 
V is not too small, although we will not attempt to provide a detailed statistical analysis on this point. 

Returning to the case of the phase-estimation procedure in Fig.[7I^ii), /d is upper bounded by a value 
termed the quantum Fisher information /qu . For the case of the interferometric detection of phase shifts 
using particle counting, /qu can be straightforwardly derived (see Ref. HOj a special case of the analysis 
in Ref. '2T'), and for the interferometer in Fig. [71[ii) the value is /qu = 4A^m7 (where A^my denotes the 
particle-number variance (wy) — (my)^). It is important to verify whether this bound is tight across 
the range of possible values for x (many existing schemes are in fact suboptimal in this respect). To 
check this, we verify that the approximate NOON states at 7 and 8 satisfy the path-symmetry condition 
identified in Ref. 20. Applying the analysis to the amplitudes for positions 7 and 8, after toi, mi and mg 
particles are counted at detectors /?i, £'2, and /)g, 

^mi ^"^2 qJ^"^ (l^^ 

Cmum2,rmms,mg OC (0| / ,^ ,^ ,^ , l^aNjs) (38) 

V !TO2 ''TlT Img iTOg ! 

the path-symmetry condition requires that, 

(39) 

where the phase factor e^'^ is the same for all possible values for 7717 and mg = N — M — my (where 
M = mi + m2 -f- mg and N ~ + A^/3), and indices 7 and 8 have been swapped. For the amplitudes 
C'mi,m2,mT,ms,mg, this Condition can be verified by inspecting the explicit forms of the operators a.^, given 
by Eq. ([2]), under (scalar) complex conjugation. As previously it is assumed that {9 = 7r/2,ie*^ = ±1 
and C ~ tt/2). Since ai, 02, and og are proportional to Oq ± a^, complex conjugation generates only a 
fixed phase factor contributing to 7 in the path-symmetry condition. In addition we find (ag)* = iar 
and (ay)* = iag. Hence complex conjugation swaps the 7 and 8 labels, and contributes a fixed factor of 
(my + mg) 7r/2 to 7. 

To compare the usefulness for phase estimation of states with the same total particle number my -l-mg, 
we adopt the quality factor of Eq. (|37|) qi — 4: (A^my) / {rriT + m^)'^-. qi is normalized between and 1, with 
the maximum value being attained by a perfect NOON state. For the input state {NaNp) combined at a 
beam splitter, with Na — Np — {rri'j -l-mg) /2 (assumed to be even), qi has the value 1/2 -I- l/(my -I- mg). 
In particular, there is already a Heisenberg-limit type scaling for the precision, which at larger particle 
numbers is less than that for a corresponding NOON state by a constant factor of \/2. 

While the direct dual-Fock-state method has a smaller precision by \/2 than the corrected NOON 
method, one must also consider how many particles the corrected device "loses" in the detectors Dl, D2, 
and D9, especially if the total number of quantum sources N is limited. If the final output contains a 
fraction / of the original source number, myg = fN, then to exceed the dual-Fock-state method in phase 
accuracy we require 

< V2 (40) 

For (j2 — 0.95 we have / > 0.72. For, say, = 60 that would require myg > 44. We have repeated the 
averaging calculation shown in Fig. 9 for myg = 44 and find in that case qi — 0.81, qi = 0.96 (still a very 
high value!). However the absolute probability of that particular value of myg is only 6 x 10~^ compared 
to the value 0.02 for myg = 20 of Fig. 9. Higher values of myg become even less likely. Thus the direct 
dual-Fock-state process may well provide a more efficient use of limited source resources. 



B. Probing the state 



Here we consider the uncorrected case of Sec. IIIII where we took mi = m2 so the NOON is formed in 
arms 5 and 6 and interferes at the last beam splitter with detection in 7 and 8. We take C = we also 
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have 9 = Ti/2 and ^ = (as above in Sec. in which case we find 



K 



mi ,m2,m7 ,ms — i i 



p!(mi - p)\{Na - p - m8)!(m2 + mg - + p)\ 



(41) 



where i^T is a normahzation factor. A plot of this probabihty versus my is shown in Fig. HI The oscillations 
are equivalent to interference fringes and are similar to those found in Ref. [13 where the phase states 
shown in Fig. [TU] were allowed to interfere. 




FIG. 8: Plot of Pmi.mamT.mg of Eq. (|4T]| versus rm for mi = m2 = 40, Na= Nfj = 40. Here C = The solid line 
is the exact result; the dotted line is the approximation of Eq. (|4ip . 



The NOON superposition in the two arms is analogous to macroscopic particle interference in a two-slit 
experiment. As explained in more detail Ref. [l^ any attempt to detect the arm in which the N particles 
travel results in the destruction of the interference pattern. 



VIII. CONCLUSION 



In conclusion, we have discussed a measurement-based approach to gene rating atomic NOON states 
with a high particle number, developing ideas first proposed in Refs. [l3. [l5l . The key requirements 
for the method are dual-Fock state Bose- Einstein condensates for the input, atom interferometry, and 
particle counting, the basic experimental feasibility of which have already been demonstrated [24, 
In the photonic case both high-efficiency number-resolved detection [26] and feed-forward switching [27| 
have been demonstrated experimentally. However photon loss would make successful implementation of 
the scheme in a truly scalable manner very difficult [28,]. We have shown using two different NOON 
quality factors that, when the number of detected particles is sufficiently high, the NOON states at the 
output are a very good approximation to the ideal. While dual-Fock states assumed for the input enable 
measurement precision better than the standard-quantum limit, and Heisenberg- limit-like scaling, NOON 
states saturate the Heisenberg limit, and represent the optimum strategy using nonclassical resources. 
One can also envisage macroscopic NOON states being used to demonstrate the quintessential two-slit 
experiment, in analogy to a large molecule propagating in a superposition state through a double slit, 
and re-interfering with itself at a screen. 

In terms of the detailed analysis, the methodology differs somewhat from that used in Ref. [ij. In 
contrast to the previous work, an integral representation of (5-functions is used (rather than working in an 
over-complete basis of coherent states). This approach is simpler; for example, the interferometer parame- 
ters can be determined by inspection of the binary form of the annihilation operators at various positions. 
We derive rigorously an algebraic condition for the transmission coefficient for the side-detections for the 
correction stage. The optimum value for the coefficient cannot be attained in practice, since it depends 
on the unknown outcome at the side detector itself. This problem is solved by substituting the mean par- 
ticle count - the best strategy possible. A simple expression for the transmission coefficient is derived by 
considering the mean intensities in different arms of the interferometer. This leads to a new explanation 
of the feed-forward method, and the achievability of high values for the NOON quality factors. 
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Appendix A 

We present here the derivation of Eqs. ([5^ and We start with Eq. (|28p and write 



1 



4"TOi!m2! 

Then for fixed m^i + mg + mg = N — mi — m2 we have 



m^f ,7nQ,7nQ 



777,5/ 'mgimg! 



(42) 



(43) 



where A/" is a normalization constant. Change summation variables from m^' to p — m^i + mg. Then we 
can pull the Gp^me factor out of the sum on mg to get 



mgip! 



E 



mg 



mg!(p - mg)! 



Because T + R — 1 the sum on mg yields simply pR to give 

{mg)=MRj2\ rrJ- r:Gp+^,iM - me) 

An exactly analogous calculation leads to 



(44) 



(45) 



(me 



Op+mg.M 

mglpl 



meGp+me E/ 



plJimgrpp-mg 

mg\{M — mg)! 



The mg sum is just (i? + T)p = 1 so that 



{me) 



-^e{ 



me 



me^iM - me)! 



G 



p+me 



The resulting relation is 



(mg) = (1 - T) (iV - mi - m2 - (me)) 



which is Eq. (P^)) . 

By particle conservation we have 



(ms + me) = N — mi — m2 



(46) 



(47) 



(48) 



(49) 



Note the prime is missing from the ms number here (ms = ms' + mg). 

We can also prove a relation between (ms) and (me). From Eq. ([7]), with appropriate values of the 
phase shifts entered, we have 



- '>Tli,m2 ,7715 ,771-6 



4^mi!m2!m5!me! J 27r J 27r 

, -1 m2+m(3 

X ~ l)(e^'^- 1) 



(e"*"^' + l)(e*'^ + l) 



(50) 
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If we change variables to A = (0 — and A = (</> + 0')/2 we find 

NJNsl r d\ r dA ,,,, , 

J^r«i,m2,m5,m6 = _ , , T |o/V / ^/ TTT COS [( A^a - ) A] 



m 



i!TO2!m5!m6!2^ 27r J 27r 



X [cosA + cosA]"^+"^ [cos A -cos A] (51) 

To find (ms) multiply this by ms and sum over and mg subject to the restriction that + mg = 
A^ — TOi — the sum is 

S = Y' , [cos A + cos A] "Mcos A - cos A] ""^ 



(cos A + COS A) ^ 



(2cosA)^-"'i-"=-^ (52) 



This results in the average 



(A^ - mi - m2 - 1)! 
(ms) = ATmimj (mi + l)F(mi + 1, m2) (53) 



where 



dA r dk 
mi!m2! 27r 27r 
x(cosA)^-"'i-™^ [cosA + cosA]™' [cos A - cos A] (54) 



F(mi,m2) = f / / — cos [(A^„ - A^^) A 



and Kjnim2 contains other factors unimportant for our purposes. A similar equation is found for {me): 

(me) = K„r^rn2 {^2 + l)F{mi , TO2 + 1) (55) 

For large A^ — toi — m2 the (cosA)^-™i-™2 factor peaks very sharply at A = and can be replaced 
by Dmim2^iA), where Dmim2 is a factor that can be lumped into Krn,im2- The result is that 



F(mi,m2) = 1 — r / ^ [1 + cos A]"' [1 - cos A]" 

mi!m2! 27r 

f 2"'i+"^+ir(mi + i)r(m2 + i) 



yielding the result 



which is Eq. (PT|) . 



mi!m2! (mi+m2)! 



(tos) ^ mi + i mi 
(me) m2 + ^ m2 



Appendix B 



(56) 



(57) 



We want to find a rigorous formula for (mg) . Starting from Eq. (|28p we expand the operators, take 
matrix elements, yielding a (5— function, which wereplace by an integral: 

0™„™„™,„™„™, - 2^V mi!m2!m5-!m6!m9! 

X r ^e-^^°* (e^* + i)™^+'"='+™« (e^0 _ i)™=+™^ (gg) 



(59) 



-\/mi!m2!m5'!me!m9! 

/ / \ mi+mcz+mg / i\m2+m& 

* c„4 .4 (60) 
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If Na = Np the integral can be done analytically. We sum the result over all 7715/ and mg to give: 

N-M 



mi ,m2 ,7719 



E 



xr 



1 + mi + ?7i5' + mg 



1 + 7712 + - - 7715' 



(61) 



where A/" is a normalization factor. 

Fig. [9] shows a plot of Pmi,m2,mg of Eq. (|6T|) for a set of variables having a large value of (mg). To get 
this plot we used T = 1712 /mi. We find the exact result (mg) = 53.6 compared to the approximation of 
54 given by Eq. (|33)) . 

Prob 



0.06 
0.05 
0.04 
0.03 
0.02 
0.01 



20 




FIG. 9: Pmi,m2,ma of Eq. gS)) versus mg for Na ^ ^ 70, mi = 40, m2 = 10. We find (mg) = 53.6, while the 
approximate formula of Eq. (|33p gives 54. 



Appendix C: Circuit Efficiency comparisons 

It is likely, in any set of experimental runs, that a random assortment of values of mi, m2, and mg will 
be averaged over in making a NOON state. What percentage of the inputs will result in good NOON 
states? Here we compare the corrected and uncorrected efficiencies to get an idea of how successful is the 
correction circuit in providing good NOON states. 



A. Uncorrected Circuit 



With the uncorrected circuit one can, for a given N value, compute in the 2D space of {mi,m2}, the 
various possible NOON output numbers, msg — N — mi — m2, and the corresponding NOON quality 
factor qi and total absolute probability (probability normalized over all five m variables) of getting each 
result. In a 2D plot of the probability one finds the highest qi factors along mi — m2, as expected. This 
line of high qi is a minimum of absolute probability - nonequal mi , m2 values are more probable. A 
sample is shown in Fig. [TU] for iV = 60 and msg = 20. Note that only the middle seven mi values give 
a quality factor greater than 0.90 and only the middle three are greater than 0.95. When msg is larger, 
the quality factors drop so that for, say, msg = 30 only the case mi = m2 = 15 gives qi = 0.95 and 
only the middle three points have qi > 0.90. As msg decreases, the distribution of high quality factors 
broadens greatly, but one is then getting fewer NOON particles as output, and the absolute probability 
of occurrence of these states is lower. 

We can average the data of Fig. [10] over all mi to get a NOON that would occur if we accepted all events 
at constant msg = 20. The result is shown in Fig. [TT] The resulting quality factor is low (0.26) because 
we have averaged over some very poor states. The total probability of getting any results corresponding 
to msg = 20 is 0.0036. 

We can assume that someone using the uncorrected circuit would be more selective and keep only data 
associated with mi values that have higher quality factors and sufficiently large m5g values. Thus if one 
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10"* Prob 




FIG. 10: Uncorrected circuit quality factor (left) and absolute probability (right) for fixed A'^ — 60, mse = 20 as 
a function of mi. The best quality occurs for mi — m2. Note that the probability is smallest for the best quality 
factor. 




FIG. 11: The NOON state resulting from averaging over all mi corresponding to fixed A'' — 60, msg = 20. The 
quality factors are low: gi = 0.27 and q2 = 0.53. 



demands that the NOON output particle number msg be 20 or larger (still for N — 60) for qi — 0.90, 
one will get this 6.2% of the time. Table |TT] shows the percentages for various output minimum numbers. 



B. Corrected Circuit 



In the corrected circuit we have one more variable mg, making it very difficult to get the data set 
equivalent to that which led to, say. Table |TT1 However, we can analyze a particular case corresponding 
to fixed values of N and output number myg. We again choose = 60 and myg = 20. First we present a 
NOON that is a probability- weighted average over all TOi,m2, and mg corresponding to those values. Of 
course, we have mg = N — myg — mi — m2 determined for each set of the variables mi and m2. Thus we 



TABLE II: Probability -Pjv^;,, (in percent) for getting a NOON with quality qi and output number greater or 
equal to Nmin for A'^ = 60 for an uncorrected circuit. 

A^m^n % (gi = 0.90) %(gi = 0.95) 

35 0.30 

30 2.7 0.2 

20 6.2 1.6 

15 6.2 2.0 
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have a double sum with the probabihty distribution given by 



Ml 2 



mi— 



Af-mi 
m2^mi + i 



m2,m7,m8,77i9 



(62) 



where M — N — mrs and the factor of 2 in the second term takes account of the symmetry that occurs 
when mi and m2 are interchanged. The NOON state produced by this process is shown in Fig. 1121 
The efficiency in getting this high quahty resuh {qi = 0.94) is good; that is, output corresponding to 

Prob 




10 



15 



20 



OT7 



FIG. 12: The NOON state resulting from averaging over all mi and m2 corresponding to fixed TV — 60, 
The qi quality factor is 0.94, while 52 = 0.98 



rriTs 



20. 



any element with this output number mys occurs with absolute probability 2.1%. This NOON quality 
comes without any selection of specially chosen values of mi and m2 as occurred in Table HI The point is 
that the mg probability distribution is small away from the high quality points because the transmission 
coefhcient has been chosen properly. 

We can see how the probability distribution selects high quality by taking apart the above NOON 
output. An important aspect of this is the position in {mi,m2} space for which the actual value of 
mg — N ~ myg — mi — m2 is equal to its most likely value, given by Eq. p3p . This is the position in 
the space where the absolute probability will be the largest. Moreover, we expect that the quality factor 
will be largest here too. Fig. [T51 shows the plot of the points where this match of the real mg agrees with 
its most probable value. Suppose we now pick a value of mi and plot the quality factor and probablity 




FIG. 13: Plot of the position in the {mi, 7712} grid at which mg matches its most probable value for A'^ = 60 and 
m78 = 20. 



versus m2: Fig. [HI shows these for mi = 10. There are two peaks in the quality factor corresponding to 
where m2 crosses the places of maximum probability seen in Figure 1131 The cusp occurs at the point 



17 



mi = 1712 = 10 where one switches from having T = mij'mx to the inverse. The probabiUty shows only 
the peak at the second crossing; there is also a peak at the first crossing, but it is too small to appear 
on the graph. Other samples of such plots show the probability similarly peaking at the point of largest 




FIG. 14: Plots of the quality factor q\ (left) and absolute probability (right) versus vn-2 at fixed mi = 10 for 
= 60, m78 = 20. 

quality factor. 

Finally we plot in Fig. [12] the quality factor and probablity, summed over all r7i2, as a function of mi, 
again for fixed N = 60, m7g = 20. There are two peaks in the probability since, according to Fig. [T51 
there are equal probability regions in the {mi,m2,} space symmetrically at {10, 20} and {20, 10}. 

*i Prob(%) 




FIG. 15: Plots of the quality factor q\ (left) and absolute probability (right), summed over all m2, versus mi at 
fixed mi = 20 for iV = 60, myg = 20. 

Clearly the correction process is successful in producing good quality NOON states with high proba- 
bility. 
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